Optical bistability in sideband output modes induced by squeezed vacuum 
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We consider N two-level atoms in a ring cavity interacting with a broadband squeezed vacuum 
centered at frequency cu 3 and an input monochromatic driving field at frequency ui. We show that, 
besides the central mode (at ui), an infinity of sideband modes are produced at the output, with 
frequencies shifted from uj by multiples of 2(lu — uj s ). We analyze the optical bistability of the two 
nearest sideband modes, red-shifted and blue-shifted. 

PACS number(s): 42.65Pc, 42.50.Dv 



I. INTRODUCTION 

Optical bistability (OB) has been the subject of intense research since its prediction and observation in the 1970's 
[1-3]. In Ref. [4] a model consisting of a system of homogeneously broadened two-level atoms driven by a coherent 
resonant field proved to give a successful description of OB. Due its potential applications in optical devices there has 
been a lot of efforts to observe and understand the phenomenon of optical bistability in two- level atoms [6-18]. 

The effects of the squeezed vacuum field on the absorptive OB for a system of two-level atoms in a ring cavity (see 
Figure 1), with different relaxation rates of the in-quadrature and in- phase components, were originally calculated 
in [8]. The authors verified that the squeezed vacuum strongly affects the OB, through the increase of the atomic 
decay time and through the introduction of a relative phase between the input pumping and squeezed vacuum fields. 
Although several aspects of squeezed vacuum effects on OB have been considered [9-12], no explicit calculations where 
done, to our knowledge, to the situation where the frequencies of the input fields, pump (uj) and broadband squeezed 
vacuum (carrier lo s ) are detuned. In papers [9-12], exact resonance between pump and squeezed fields frequencies, 
u; = lo s , were assumed in order to maximize the squeezing effects. Nonetheless, consideration of detuning, u> ^ u> s , is 
the source of interesting physics as to be shown in this paper. Here, we analyze the effects of that detuning over the 
OB in the output field, produced by a system of two-level atoms in a cavity. We show that, besides the central mode 
at u), the output field contains an infinity of sideband modes at frequencies shifted from to by multiples of 2(w — oj s ). 
We analyze the OB of the two nearest sideband modes, red-shifted and blue-shifted. 

The paper is organized as follows: In Sec. II we introduce the model we use, and derive the system dynamical 
equations. In Sec. Ill we obtain the stationary solutions for the output field. In Sec. IV we discuss the results and 
present our conclusions. Finally, in Appendix A we derive the many-body master equation and apply the mean-field 
approximation for a dilute atomic gas. 

II. MODEL 

We consider an input pump coherent signal of undepleted electric field amplitude Ei n and a broadband squeezed 
vacuum, with frequency distribution centered at oj s , interacting with TV two-level atoms. The Hamiltonian of the 
system is given by 

H = ~w S + F*e iut S- + Fe-^S+ + 5> fc #& fc + ]T (g k b k S + + g* k b+S-) , (1) 

k k 
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(we have set h — 1), where, the first term stands for the two-level atomic system (transition frequency u>q), the two 
following terms represent the interaction between the atoms and the input pump field amplitude Ei n , F = jiEi n ( p 
is the atomic dipole moment), the fourth term corresponds to the squeezed vacuum modes and the last one is for the 
interaction between atoms and squeezed vacuum field. Operator bk (b^) annihilates (creates) squeezed field quanta 
of frequency LOk and is the coupling constant. The atomic collective operators are 



N 



(2) 



i=l 



where so(i) and s±(i) are single particle operators satisfying the commutation relations [so(i), s±(J)] — i26%js±{i) 
and [s+(i), «-(.?)] = Si,jSo(i). Although the atoms do not interact directly with each other and the coherent field is 
assumed undepleted, they become correlated to each other, only due to their coupling with the squeezed vacuum field. 

In the mean field approximation and in a rotating frame at frequency u), the atomic system is described by an one- 
body master equation, obtained by calculating the trace over the squeezed vacuum degrees of freedom (see Appendix 
A for a detailed derivation), 



dp(t) 
dt 



= - [H eff ,p(t)] - | e l9 c let Q - iv} 2 sinh r cosh r s+ps + + h.c. 

+ ^ sinh 2 r (s-S + p — 2s + ps- + ps_s + ) + ^ cosh 2 r (s + S-p — 2s_ps + + ps + s_)| . 



(3) 



The term in braces represents the phase-sensitive damping due the squeezed vacuum, r is the squeezing parameter, 
e = 2 (lj — u> s ) is twice the detuning between input pump and squeezed vacuum fields, 7 is the damping constant 
and 9 is a phase reference of the squeezed vacuum field. H e ff is an effective nonlinear mean-field single particle 
Hamiltonian, describing the motion of one atom in the sample, 



H eff = ^ (5 - v cosh 2r) s 



N - 1 



+ h.c. 



(4) 



where (s±) = Tr(ps±), S = uo — u is the detuning between atomic transition and pump field frequencies and v is the 
dynamic frequency shift; being much smaller than 7, it will be neglected [19]. From second term in the hamiltonian 
(4) we see that effectively a single generic atom is excited by the imput field plus a polarization field 



N — 1 / 7\ 
e po l(t) = — — {-" + i 2) ( s +) 



due the other (N — 1) atoms. 

The equations of motion for the atomic operators mean values are 



(5) 



(s ) = 2i/i (e T (t)<s-> - 4(*)<s_)*) - 7 ((so) cosh2r + 1) , 



(6) 



(si) = -ifi(*_) + ipe* T (t)(s ) - Qe 4et (s_)*, (7) 

where (s+) = (si)*, n = S - i( 7 /2) cosh2r , Q = (~f/2)e t0 sinh 2r, and 

e r (t) =E in + e po i(t) (8) 

is the total effective field experienced by a single atom. The second term in f2 is due to the commutation relations in 
the Heisenberg equations. Furthermore, in the induced atomic polarization field 

e poJ (t)^A(s_)*(t)/ M , (A = iJjV e// ) (9) 

we have assumed an effective number of atoms N e ff contributing effectively to this field ( N e f f <c N) . 

In the next section we obtain the stationary effective field amplitude er{t) as function of E^ n and system parameters. 
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III. STATIONARY SOLUTIONS 



For no detuning between vacuum squeezed and pump fields, e = 0, there is no explicit time dependence in Eq. (7), 
and the equilibrium solutions ((s'_) = 0, (s'o) = 0) (s_) e<? and {so) eq are easily obtained as function of the output field 
ct, which, together with Eq. (8) enables to recover the well known result [8,10], 

E in = e T : 7M (^-y4) . (10) 



Vfi/ |e T | 2 - 2^ 2 (Qe 2 T + Q* (e^) 2 ) - 7 (|ft| 2 - |Q| 2 ) cosh2r ' 



The bistable behavior becomes evident from plotted output field amplitude modulus \ex\ as function of the same for 
the pump field \Ei n \, as displayed in Figure 2. Above a critical value of N e ff an 5-shaped curve is produced, meaning 
that there are two possible output fields for a single input one. Moreover, the S*-shaped curve is quite sensible to the 
phase difference between input and squeezed vacuum fields, as stressed in Refs. [8,10]. 

We are interested in the situation e ^ 0, when equations (6)- (7) arc no more autonomous, so the asymptotic 
stationary solutions are periodical time-dependent series 

oo oo 

( S _}= vj a « emet ' ( s °)= E 6 « emet ' ( n ) 

n— — oo n— — oo 

whose coefficients can be determined from equation (6) and (7). The output field amplitude is also expanded as an 
infinite series 

oo 

er(*)= E £ ^ ntt - ( 12 ) 

n— — oo 

In a non-rotating frame, the total output field amplitude is a superposition of an infinite and countable number of 
modes, 

E T {t) = e T (t)e- wt = £ e- wt + £ +1 e~ 1 ^'^ + £_ 1 e- i <"' +e > t + ... , (13) 

at frequencies u) n = u> ± ne, for n = 0, 1, 2, .... 

Inserting the series (11) and (12) into Eqs. (6)-(7) and equalling coefficients with same time dependent factor e me *, 
one gets the following equations in terms a n and b n , 

E ln =£ Q - Aa* , for n = 0, (14) 
£ n = Aa*_ n , for n + 0. (15) 



% (ne + Q) a n + Qa*_ n+1 = i £* n (16) 



(me + 7COsh2r)6 n 

oo 

= 2i E ( £ n- m a m - £* n - n a*_ m ) - j5 nfi - (17) 

m— — oo 

After a lengthy but straightforward algebraic manipulation of equations (15)-(17), one obtains an equation involving 
only the coefficients a n and the central output field amplitude £q, 

£* ( £* 

G n {e)a n + F n (e)a*_ n + Qa*_ n+1 +i —2— 5 nfi = -2 [ yj^ [Aa*_„aj - A*a n -ia*_ t ] 

^* ^ p i I 

+ tptt [£ a n -ia n - £*a n -ia*_ n ] + ^ v , \ \^\ a n - m ai- n ai - (A*) 2 a n - m a n -ia*_ l > (18) 
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where 



G n (e) 
F n (e) = 



ne + fl + 

, (gp) 2 
X(e) 



Yo(e) 



, \£oV 



and 



Y n (e) = ine + 7 cosh 2r. 



(19) 
(20) 

(21) 



On the left-hand-side (LHS) of Eq. (18) N e ff enters only in G„(e), while on right-hand-side (RHS) it enters the terms 
involving the products of a„'s. For field intensities of sideband modes much weaker than the central mode, we neglect 
the nonlinear terms on the RHS of Eq. (18). This allows us to rewrite the LHS in terms of a finite difference equation 
for a„, 



B n (e)a n + C„(e)a„+i 
E (e)S nfl + iJi(e)<5„ 4 , 



where 



B n (e) 
C n (e) 
D n (e) 
E n (e) 
H n (e) 



G„(e) 



F n (e)F* n (e) 



D n (e)a n ^ 



G*_ n (e) G*_ n+1 (eY 



Q*F n (e) 
QFt n+1 {e) 

G*- n +l( € ) 



cosh 2r 



£oF n (e) 
G*. n {e) 



G*_„ +1 (e)cosh2r' 



(22) 

(23) 
(24) 
(25) 
(26) 
(27) 



Even in this very linear approximation the n-dependence in the coefficients (23)-(27) does not allow obtaining an 
exact closed solution to Eq. (22), for e / 0. In the present analysis, we are going to determine only the first three 
sidebands coefficients ao and a±i. From Eq. (22) one gets the following system of equations 



B a + C ai + £> a-i = E 
Bxai + Cia 2 + Dxao = Hi 
B_ia_i + C_ia + £>_ia_ 2 = 0, 



(28) 



which is not closed because ao and a±\ are coupled to a±2, that, by their turn, are coupled to higher order coefficients. 
Instead of simply disregarding a 2 and a_ 2 in Eqs. (28), we consider a better approximation by estimating them from 
truncated continued fractions. Setting 



y- 



-(«+!) 



CL n — \ CI n 

for a n _i 7^ 0, a_ n ^ and n ^ 0, 1 we can write Eq. (22) as two equations, 



-Dn 

B n + C n x n+ i 



n ^0,1, 



(29) 



(30) 



»-» = b — ' (31) 

for positive integers n. For n = 2 in (30), n = 1 in (31) and truncation of the continued fractions, up to a second 
order iteration, yields (a higher order iteration does not affect significantly the result) 
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Substituting ci2 = x 2 cii, a-2 = U-id-i, we get the coefficients 




Bi+Ci4' 



(34) 



ax = 



Hi - Pi a 

B 1 + Ci4 2) 



a_i = 



(35) 



IV. RESULTS AND CONCLUSIONS 



Using the solutions for the amplitudes oo and a±i, Eqs. (34) and (35), we can analyze the functional dependence 
of the output fields £ and £± as function of the input field E in , in modulus. To simplify the illustration of the 
bistable behavior, we assume 9 being the phase difference between pump and squeezed input fields. The output field 
amplitudes |£o|, = |^ a -i| an d = |Ao*| are plotted as functions of E in in Figs. 3-(a),3-(b), and 3-(c) 

respectively. The parameters are set as N e ff = 101, e/7 = 2.0, r — 0.5, 5 = 0, and 9 = tt. We verified that the OB 
looses the phase-sensitivity, varying more significantly with r, because the coefficients a , a\ and a_i now depend on 
\Q\ 2 , instead on Q. Although the sideband field intensities are much weaker than the central one, they also display 
a bistable behavior, with turning points occurring at the same input field intensity. The dashed lines correspond to 
the unstable branches, the arrows indicate the path followed by the output field variation as the input is increased 
or decreased. The bistable behavior of the central mode (Fig. 3-(a)) is similar the case where e = 0, however, the 
sideband modes, £\ (Fig. 3-(b)) and £-\ (Fig. 3-(c)), that are respectively, red-shifted and blue-shifted with respect 
to the central mode, show some qualitative differences. Differently from the central mode, at strong pump amplitude 
modulus, the sideband modes show a monotonic decrease in the amplitude modulus at the output. The sideband 
modes also present the following different features in the switchings, or jumps from low to high amplitude (and vice- 
versa) in comparison with the central mode: i) By increasing the input field intensity the (a) — > (b) switch is from 
low to high amplitude, in modes £q and £+1, see Figs. 3-(a) and 3-(b), however it is inverted in mode £-1, switching 
from high to low amplitude, see Fig. 3-(c). ii) By reverting the path, going from high to low input intensity the 
switches occur from high to low output amplitudes, (c) — > (d), in modes £0 and £ + \, Figs. 3- (a) and 3-(b), while it is 
from low to high in mode £-\. Essentially, the sideband modes show inverse behavior, with respect to the switchings, 
iii) Comparatively to the central mode, the sidebands present a higher contrast in the jumps from higher to lower 
amplitude. 

A possible application of the above results could be the simultaneous transmission of a message by the output 
field through three different channels (the three modes), where the triplicated information could be useful for error 
control. Additionally, the codification in the blue-shifted sideband (0,1,1,0,0,1,...) is the inverse of that in the other 
mode (1,0,0,1,1,0,...), so the sideband modes could transmit information as like the codification occurring in the DNA 
double-strand macromolecule, where one strand sequence is the inverse of the other. 

In conclusion, we have shown that the frequency detuning between input pump and squeezed vacuum fields, inter- 
acting with two-level atoms, gives rise to a multiple-mode output field with frequencies that are multiples of e. By 
analyzing the closest (red-shifted and blue-shifted) sideband modes, to the central one, we did verify new features in 
the bistable behavior. Although the obtention of a pump and squeezed fields with controllable phase difference could 
be, at the moment, experimentally difficult, because both should derive from a common source, we believe that the 
reported physical effects could be useful in optical devices and in the transmission of information. 
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APPENDIX A: MASTER EQUATION 

In a referential frame rotating at frequencies u — ujk the hamiltonian (1) becomes 

H = H 0S + V(t), (Al) 

where 

H QS = ^So + F*S-+FS + , (A2) 

the collective operators are defined in (2), S = u> — ujq and 

V(t) = ]T (gkbkS+e*"-""* + h.c) . (A3) 
k 

Following the usual procedure [19], by eliminating the reservoir degrees of freedom one obtains a pre-master equation 
for the system density operator p(t) 

^ = \ [H 0S ,p(t)] J* dt'Tm \V(t), [V(t'), p(t')p R }} (A4) 
where pn is the state of the reservoir, at thermal equilibrium. Substituting the interaction (A3) in (A4) one gets 

^ = -i [Has, p{t)]~ f ' dt' {£n(t,t') [S+, [S+, p(t')]] + £ 2 (t, f) [S-,S + p(t')} (A5) 
H2i(t,t')[S + ,S-p(t')] + h.c.}. (A6) 
The coefficients £ij(t, t') are characterized by the kind of reservoir, 

f dt'£n(t, t'W) = [ dt' V g u ,g u „e i ( u - u ') t+i ( u - u ") t ' (b u ,b^) R p(t% (A7) 
J o Jo u/u „ 

f dt'£ 12 (t,t'W) = f dt' g^e-<—') t+ <—"y (b+b LU „) R p(t / ), (A8) 

f dt'^t'W) = fdt 1 J2 9»>9l»e<"-" ')*-<— ")*' (b u ,b+,) R ptf). (A9) 

where r = Tr^ (pB.b w ib w "). For a squeezed reservoir 

(b u 'b u ") R = -e i9 sinhrcoshr 8[u" - (2w a - w')] 5 (A10) 

(b^b u n) R = sxnh 2 r 5{w'-J'); (b u >b+„) R = cosh 2 r S [J - w") (All) 

where r is the squeeze parameter, 6 is the reference phase of the squeezed field, and uj s is the central resonant frequency 
of the squeezing device. Going from sums to integrals in Eqs. (A7-A9) and using expressions (AlO-All), one gets for 
example 

/ dt'£u(t,t')p(t') = -e^sinhrcoshr / dcj'D(u)')g(u')g(2u) s - u')e i( - u - u '^ / dt'e^' 2 "' p(t') 
Jo Jo Jo 

where D{u>) is the reservoir density of modes. Making the change t — t' = r and invoking the Markov approximation 
p(t — t) ~ p(t) we obtain 



I d Te -H<>>-z*.+v')T - p q. - T )S£ dre-*"-*"'*"")* p(t) = p{t) 
Jo Jo 



nS (u) - 2uj s + uj') - iV 



1 



uj — 2uj r + oj' 



and 



with 



/ dt'Zn(t,t')p(t') S| n (t)p(t) 
Jo 

r r°° 

| n (t) = -e 4e sinhrcoshr e 2i ("-^)* n D(2u a - uj)g{2uj s - oj)g{uj) - VP \ duo 

Jo 



,D( uj')g(Lj')g(2Lj s -u') 
uj — 2uj s + uj' 

where V stands for the Cauchy principal value. For \uj s — u>\ <C uj the two terms in the brackets are assumed being 
approximately constant, so we define the damping constant (7) and the dynamical frequency shift (V s ) 



7 = 2nDg 2 , 



v = V [ 
Jo 



doj 



,D{uj')g{uj')g{2uj s -uj') 



oj — 2uj s + uj' 
therefore 

fn(t) = -e 4e sinhrcoshr e 2i (<"-^)* (J_ _ ^ 
Following the same procedure one obtains the other coefficients, 

I12 = Q - w s j sinh 2 r; £21 = Q - if) cosh 2 r, 

which are time- independent. 

Thus the master equation for an iV-atom system becomes 

= \ [H$\p N (tj\ {| n (t) [S+, [S+, p N (t)]] + & S +Pjv (t)] + 61 S-Piv(t)] + ft.c.} 
while for a system of p— atom system, p < N, it is 



dt 



H$,P P (t) ~ { fn(t) E [s + (i),[s + (j),P P (t)}} 



(A12) 



+ 61 ( J2 [*+(*). «-C/)Pp(*)] + - p)Et s +W' Tr P+i S-(P+l)P P +i(i)] ] +A.c 



i=l 



(A13) 



For a dilute system the atomic correlations may be disregarded, so, we shall consider a single generic atom (p = 1) 
moving in a mean field produced by all the others, with the 2— atom density operator factorized as pi « p\ <E> p\. In 
this approximation equation (A13) reduces to 



d P i{t) = 1 
dt i 



H$,Pi(t)] - {|n(t) [s+, [a+,pi(t)]] + & ([*-,*+Pi(t)] + (JV- 1) [a_,pi(t)]) 



+ | 21 ([«+, + (TV - 1) <«_) [«+,pi(t)]) + /i.e.} . 

with the single particle Hamiltonian 



(A14) 
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o-(l) 
n OS 



-s + + Fs+ 



(A15) 



and (s±) = Tr(s±pi) is the mean value. Rearranging the terms in Eq. (A14) and dropping the subscript 1 in pi we 
can write Eq. (A14) as 



dp(t) 
dt 



= -i [H eff ,p(t)} - {2 [ e * e 2i(«-«.)t Q - ii^j sinhr coshr s+ps+ + h.c. 



+ ^ sinh 2 r (s_s + p — 2s + ps_ + ps_s + ) + ^ cosh 2 r (s + s_p — 2s_ps + + ps + s_) j . 



(A16) 



The single particle effective Hamiltonian in Eq. (A16) is given by 

H eff = ^(S - j/cosh2r) s + {F*8- + Fs+) + (N - 1) [(-1/ + i^j (s+) s_ + /i.e. 



(A17) 



it contains nonlinear terms due the mean-field effect of the remaining N — 1 atoms. The frequency shift ^cosh2r is 
due to the interaction with the reservoir. The second term in the RHS of Eq. (A16) stands for the dissipative part 
due to the decay in the squeezed vacuum [19]. By setting u) — uj s , identifying sinhr coshr — ► m, sinh 2 r — ► n and 
cosh 2 r — ► n + 1, the dissipative term of the master equation takes the same form as considered in [8]. 
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FIGURE CAPTIONS 



Figure 1. Schematic ring cavity with N two- level atoms in a cell. Input signal at the left, output at the right and 
injection of squeezed vacuum from above. Mi to M4 specify the mirrors. 

Figure 2. Output versus input field amplitudes. Sensitivity to the phase 9 is manifest. 

Figure 3. Output modes versus input field Ei n . Dashed lines are for the unstable branches. Arrows indicate the 
direction of variation of output amplitudes with increasing (decreasing) Ei„. (a) central mode amplitude Sq. The 
jumps goes from (a) — > (b) ( (c) — ► (d) ) increasing (decreasing) the output amplitude, (b) red-shifted sideband 
mode The jumps are in same direction as in Figure (a), (c) blue-shifted sideband mode £-\. The jump goes 

from (a) — > (6) ((c) — > (d)) with decreasing (increasing) amplitude of the output field. Both jumps occur in direction 
opposite to those in Figures (a) and (b).The parameters are set as N e ff = 101, e/7 = 2.0, r = 0.5, (5 = 0, and 9 = it 



9 



